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Abstract. The 15 N(p,7) 16 reaction controls the passage of nucleosynthetic material 
from the first to the second carbon-nitrogen-oxygen (CNO) cycle. A direct 
measurement of the total 15 N(p,7) 16 cross section at energies corresponding to 
hydrogen burning in novae is presented here. Data have been taken at 90 - 230 keV 
center-of-mass energy using a windowless gas target filled with nitrogen of natural 
isotopic composition and a bismuth germanate summing detector. The cross section 
is found to be a factor two lower than previously believed. 
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Figure 1. Energy levels of 16 relevant to the 15 N(p,7) 16 reaction at low energy, in 
keV [1] . Primary (dotted) and secondary 7-ray transitions are also shown. 



1. Introduction 

The 15 N(p, 7 ) 16 reaction (Q-value Q = 12. 127 Me V) links the CN cycle [1, 2J to the 
CNO bi-cycle and all further CNO cycles [3]. The 15 N(p,7) 16 cross section a(E) (E 
denotes the center of mass energy in keV, E p the proton beam energy in the laboratory 
system) can be parameterized [3] by the astrophysical S-factor S(E) defined as 

S(E) = <x(E)E exp(212. 85 j^E). (1) 

At astrophysically relevant energies E < 1 MeV, the 15 N(p,7) 16 excitation function is 
influenced by two resonances at E p = 335 and 1028 keV (E x = 12440 and 13090 keV, 
figure [T]), with respective widths of T p = 91 and 130 keV, both decaying predominantly 
into the ground state of 16 [4J. For the E x = 12440 (13090) keV level, 1.2% (0.58%) 
decay branching to the + first excited state of 16 at 6.049 MeV has been reported 
[I]. In addition, for the 13090 keV level, there is 3.1% decay branching to the 1~ third 
excited state at 7.117 MeV [I]- No other decays to 16 excited states are known for E p 
< 1028 keV |I] . 

The non-resonant cross section has been studied in previous experiments using Nal 
[5] and Ge(Li) [B] detectors, reporting cross section data for E p = 150 - 2500 keV [6]. 
Citing discordant normalizations between those two studies [5j E], only the data from 
one of these studies [B] have been used in reaction rate compilations [7J [S]. 

Recently, the asymptotic normalization coefficients (ANC's) for direct capture to 
the ground and several excited states in 16 have been measured [9]. It was found that 
the low-energy non-resonant yield is dominated by ground state capture [9], but the 
new ANC leads to a much lower direct capture cross section (sum of direct capture to 
all states in 16 0) than previously |6J. The new ANC values have then been used in 
an R- matrix fit [9] including also the cross section data from refs. [HI E], suggesting a 
factor two lower astrophysical S-factor than previously believed [HI 0, E] • Another recent 
R-matrix analysis concentrating on ground state capture was based again on the direct 
data from refs. [5l E], and it also indicates a much lower S-factor |10j . In view of the 
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conflicting data [3J [6] and the recent extrapolations P, [TO], new experimental data is 
clearly called for. 

The aim of the present work is to experimentally determine the 15 N(p,7) 16 cross 
section directly at energies corresponding to hydrogen burning in novae. The relevant 
temperatures in novae [HJ [12] are Tq = 200 - 400 (Tq denoting the central temperature of 
a star in units of 10 6 K), corresponding to Gamow energies [3] of -Ecamow = 150-240 keV. 
In order to obtain the new cross section data, spectra from a radiative proton capture 
experiment at LUNA that has been performed using nitrogen gas of natural isotopic 
composition (99.6% 14 N, 0.4% 15 N) have now been analyzed regarding the 15 N(p,7) 16 
reaction. 

2. Experiment 

The experiment has been performed at the Laboratory for Underground Nuclear 
Astrophysics (LUNA) in Italy's Gran Sasso underground laboratory (LNGS). The 
LUNA facility has been designed for measuring low nuclear cross sections for 
astrophysical purposes [131 HH HS1 HH H3 HH1 US] , benefiting from its ultra-low laboratory 
7-ray background [20l El] . 

2.1. Target 

A windowless, differentially pumped gas target cell filled with 1 mbar nitrogen gas of 
natural isotopic composition (0.366% 15 N [22]) has been irradiated with E p = 100 - 
250 keV H + beam from the 400 kV LUNA2 accelerator [23J. The emitted 7-rays have 
been detected in a 47r BGO summing crystal [2lj . The calorimetric beam intensity 
values are known with 1.0% precision [24J. 

The natural isotopic composition of the target gas enabled parallel experiments 
on 14 N(p, 7 ) 15 [US [25] and 15 N(p, 7 ) 16 (present work). The 14 N(p, 7 ) 15 analysis 
is already published including full experimental details [TOl [53]; the present work 
concentrates on aspects pertinent to obtaining the 15 N(p,7) 16 cross section. 

During the experiment, nitrogen gas of natural isotopic composition and 99.9995% 
chemical purity was flowing through the windowless target cell with a flux of 2 
liters/second. No recirculation was used, so the gas was discarded after one passage 
through the target. The effective 15 N target density for the present work has been 
obtained scaling the known target density (3.2% uncertainty including the beam heating 
correction [23]) with the standard isotopic composition [22]. A recent survey has found 
that >99% of nitrogen-bearing materials have isotopic abundances within 2.0% of the 
standard value [22J, which is defined to be that of atmospheric air. The 15 N content 
of atmospheric air on different continents has been found to be constant to 2.6% [26J, 
and commercial tank gas even falls within 1.0% of the standard [27]. In order to verify 
whether these findings also apply to the presently used tank gas, gas samples of the 
type of nitrogen used here and from the same supplier have been sent to three different 
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laboratories for isotopic analysis. The isotopic ratio was found to be within 3% of the 
standard. As relative uncertainty for the isotopic ratio, 3% is therefore adopted. 

2.2. 'y-ray detection 

The 7-ray detection efficiency of the BGO detector [21] has been obtained by a dedicated 
simulation with GEANT4 [28] . The simulation has been validated at low 7-ray energy by 
measurements with calibrated 7-ray sources and at E 1 ~ 7 MeV by a detailed comparison 
with the results from the previous [25J GEANT3 simulation. An uncertainty of 3.0% is 
quoted here for the probability of detecting isotropically emitted 12 MeV 7-rays. 

The GEANT4 summing detector efficiency depends, however, also on inputs from 
experiment, such as the decay scheme and the angular distribution of the emitted 7- 
radiation. If the capture does not proceed directly to the ground state, but to some 
excited state, several 7-rays may be emitted, leading to lower detection efficiency when 
compared to ground state capture. 

In order to understand the decay scheme, germanium spectra taken at E p = 400 keV 
(slightly above the E p = 335 keV resonance) bombarding solid Ti nat N targets with proton 
beam [TH] have been reanalyzed here. Experimental upper limits of 1.9% (1.8%) for the 
primary 7-rays for capture to the excited states at 7.117 (6.049) MeV in 16 have been 
derived. In addition, from a reanalysis of germanium spectra [20J taken with the present 
gas target setup at E p = 200 keV, an upper limit of 6% for the 7-ray from the decay 
of the 7. 117 MeV state is deduced. These findings are consistent with the previous 
conclusion that for E p < 400 keV, the reaction proceeds to >95% by capture to the 
ground state in 16 O [6]. 

The GEANT4 simulation shows that the summing peak detection efficiency for 
7-rays decaying through the 1~ level at 7.117 MeV is 27% lower than for ground state 
capture. The + level at 6.049 MeV does not decay by 7-emission, so capture to this 
level cannot be detected in the 12 MeV summing peak at all. Scaling these effects 
with the above mentioned experimental upper limits for the capture probability to the 
corresponding level, 1.9% systematic uncertainty for the total cross section is obtained 
due to possible capture to excited states. 

The angular distribution has previously been found to be isotropic at the E p — 
1028 keV resonance pi), and for the present analysis, isotropy has been assumed. The 
simulation shows that due to the large solid angle covered by the BGO, the detection 
efficiency is enhanced by only 4% when assuming a complete sin 2 ?? shape instead. In 
order to account for this effect, 4% is adopted as systematic uncertainty. 

2.3. Analysis of the 'j-ray spectra 

During the experiment, 7-ray spectra were taken at twelve different incident energies 
between E p = 100-250 keV. For each beam energy, two in-beam spectra were recorded: 
one with 1 mbar nitrogen gas (natural isotopic composition) in the target, and one with 
1 mbar helium gas (chemical purity 99.9999%) to monitor ion beam induced background. 




Figure 2. 7-ray spectrum recorded at E p — 150 keV. Solid red (dotted blue) line: 
Nitrogen gas in the target (helium gas, rescaled to match the nitrogen spectrum in the 
14.4 - 18.0 MeV region). Dashed green line, laboratory background, rescaled for equal 
livetime. See text for details. 



In addition, a spectrum with 1 mbar argon gas in the target has been recorded at E p = 
216 keV. Laboratory background spectra were taken during accelerator downtimes. 

The in-beam spectra can be classified in two groups, low beam energies E p = 100- 
150 keV (example, figure [2]), and high beam energies E p = 190-250 keV (example, figure 
E]). Salient features of the spectra are discussed in the following. 

At low 7-ray energies {E 1 < 4 MeV), the in-beam 7-ray spectra are dominated 
by the laboratory background and resultant pile-up. For 4 MeV < E 1 < 8.5 MeV, the 
following in-beam 7-lines are evident [20J: 

• the 4.4 MeV 7-ray from the decay of the first excited state of 12 C populated both in 
the n B(p,7) 12 C and in the 15 N(p,a7) 12 C reactions (well visible in all the nitrogen 
spectra, visible in some of the helium spectra), 

• the ~5.5MeV peak from the 2 H(p,7) 3 He reaction (visible only for E p < 150 keV in 
both the nitrogen and helium spectra), 

• the 6.1 MeV 7-ray from the decay of the second excited state of 16 populated in 
the 19 F(p,a;7) 16 reaction (visible only for E p > 180 keV in the helium spectra), 

• the 6.2 MeV and 6.8 MeV secondary 7-rays and the ~7.5MeV summing peak from 
the 14 N(p,7) 15 reaction (well visible in all the nitrogen spectra, not visible in the 
helium spectra), 




Figure 3. 7-ray spectrum, E p = 220 keV. Red solid (dot-dashed) line: Experimental, 
nitrogen gas (simulated, assuming only the 15 N(p,7) 16 reaction). Blue dashed 
(dotted) line: Experimental, helium gas, rescaled as in figure [2] (simulated with Sim. B, 
assuming only the 11 B(p,7) 12 C reaction). 



• the ~7.7MeV peak from the 13 C(p,7) 14 N reaction (well visible in the helium 

spectra, covered by the 14 N(p,7) 15 lines in the nitrogen spectra), and 

• the 8.1 MeV summing peak from the 18 0(p,7) 19 F reaction (visible only in a few 
helium spectra). 

At E 1 > 8.5 MeV, the laboratory background [20] is negligible for the purposes of 
the present study (figure [2]) • At these high 7-ray energies, the spectrum is determined 
by only two reactions: 

(i) First, the full energy peak of the 15 N(p,7) 16 reaction to be studied, visible in 
the nitrogen spectra at E 1 = Q + E m 12.3 MeV. Because of the rather smeared 
out response function of the BGO detector to high-energy monoenergetic 7-rays, a 
region of interest (ROI) from 9.7-13.5 MeV (shaded in figs. [21 E]) has been adopted. 
The probability that a 12 MeV 7-ray emitted isotropically at the center of the 
detector leads to a count in this ROI is found to be 77% in the simulation. 

(ii) Second, two peaks from the 11 B(p,7) 12 C beam-induced background reaction (Q = 
15.957 MeV), visible in both the nitrogen and helium spectra: a summing peak at 

= Q + E ~ 16 MeV and the primary [E 1 ~ 12 MeV) 7-ray from capture to the 
4.439 MeV first excited state in 12 C. (The decay of that state has been discussed 
above.) 
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2.4- Subtraction of the n B(p,7) 12 C background 

In order to obtain the 15 N(p,7) 16 cross section, the background in the 9.7-13.5 MeV 
ROI induced by the n B(p,7) 12 C reaction must be reliably determined and subtracted. 

The n B counting rate can be monitored by the yield in the 14.4-18.0 MeV region, 
where no other beam-induced lines are present. This rate varied strongly from run 
to run, also at the same beam energy, so it was necessary to derive a background 
subtraction procedure based on data in the same experimental spectrum used also for 
the cross section determination. Assuming that the place of origin of the n B 7-rays 
is the collimator at the entrance of the target cell [20], which is hit by the beam halo 
(0.5-5% of the beam current on target), the mentioned variation of the n B counting 
rate can be explained with differences in the details of the proton beam focusing from 
run to run. 

However, even for different absolute n B counting rates, the ratio between the 
«12 MeV and ~16 MeV 11 B-induced counting rates depends only on the beam energy 
(due to energy-dependent branching ratios and angular distributions) and not on the 

jBoron. 
L 12/16 



focusing. This leads to the definition of the ratio i?? OI " n 



oBoron 4 Counts (9.7 — 13.5 MeV) 
12/16 Counts (14.4 -18.0 MeV) 



At each beam energy, the quantity -Rf 2 °/i6 nas been determined experimentally from 
a monitor run with helium gas in the target (tabled]). As a check on the reliability of 
using helium as monitor gas, at E p = 216 keV, -Rf 2 /i6 has been determined with argon 
gas instead of helium, with consistent results (tablet!]). 

The experimental R^/w vames are then compared with the results of two GEANT4 
simulations called Sim. A and Sim. B. In both Sim. A and Sim. B, the known branching 
ratios and angular distribution of the 11 B(p,7) 12 C reaction from ref. [29] are included. 

Sim. A The point of origin of the n B 7-rays was assumed not to be the final collimator, 
but the beamstop (tablet!]). 

Sim. B The point of origin of the 11 B 7-rays was assumed to be the final collimator as 
discussed above (table [U figure [3]). 

For all data points, Sim. B is closer to the experimental data than Sim. A. However, 
at the lowest and highest proton beam energies the experimental -Rf 2 °/i6 vames tend to 
be even higher than the simulated ones from Sim. B (tablet!]). In order to understand 
this phenomenon, it should be noted that the simulation results depend strongly on 
the assumed branching ratios, angular distributions, and angular correlations. The 
branching ratio is known experimentally also for off- resonant energies [29J. However, 
the angular distribution is only known at the E p = 163 keV resonance [29]. It seems 
plausible that given this limited input data, the simulation does a better job close to 
E p = 163 keV than far away, at the lowest and highest proton beam energies. 

For the actual data analysis, the experimental -Rf 2 °/i6 vames have been used. In 
order to err on the side of caution and quote a conservative uncertainty on the adopted 



Direct measurement of the N(p,7) O total cross section 



CO 
-Q 

> 

CD 
CO 



500 



400 



300 



200 



100 



A 




100 



— I — 

200 
E[keV] 



300 



— r> 
400 



Figure 4. 15 N(p,7) 16 astrophysical S-factor. Experimental data from ref. [5] 
(blue circles, limited to E > 210 keV), ref. [5] (green triangles) and the present 
work (red filled squares). Error bars reflect statistical and systematic uncertainties 
summed in quadrature. Dotted line, previous low-energy extrapolation by the NACRE 
compilation [S]. Dashed line, previous R-matrix fit, and shaded area, its quoted 17% 
uncertainty [5]. 



R^2/i6 value, for Ai?f 2 ° r °g either the statistical uncertainty or ±1.0 (an upper limit on 
the full difference between Sim. A and Sim. B) was used, whichever is greater (table [T]). 

Finally, the n B background to be subtracted in the 9.7-13.5 MeV ROI of the 
nitrogen spectrum is then obtained by multiplying the counts in the 14.4-18.0 MeV 
monitoring region in the same nitrogen spectrum with the experimental -Rf 2 /i6 varue 
from the corresponding helium run (table [T|). The uncertainty due to the boron 
background subtraction has 1.8-43% effect on the S-factor data, and it dominates the 
uncertainty for most data points. Two types of runs have been excluded from the present 
analysis: Runs that show more n B background than 15 N yield in the ROI, and runs for 
which no helium monitor run has been performed. 

2.5. Further experimental details 

The effective interaction energy has been calculated assuming a constant astrophysical 
S-factor [3] over the typically lOkeV thick target, leading to 0.7-3.2% systematic 
uncertainty including also the accelerator energy calibration [23J uncertainty. All 
systematic uncertainties are summarized in table El 
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3. Results 

Based on the spectrum integration discussed in the previous section, the 15 N(p,7) 16 
cross section has been determined at twelve effective center-of-mass interaction energies 
E eS between 90 and 230 keV (table [3]). The statistical uncertainty is typically well below 
10%. 

The present S-factor data (figure HJ) are about a factor two lower than the previous 
data by ref. [6], but still consistent at 2a level given the previous high uncertainties. In 
the limited overlapping energy region, the present data seem to agree with ref. [S], if 
ref. j5]'s low-energy data points (affected by beam-induced background) are excluded. 
The data from the present work extend to energies lower than ever measured before 
and are significantly lower than the low-energy extrapolation adopted in the NACRE 
[8] compilation. 

The present data are on average 20% lower than, but given the previous uncertainty 
still consistent with, the recent R-matrix fit based on an ANC measurement [9]. They 
are also lower than the fits shown in ref. [TU] . These R-matrix fits [21 [TUJ had relied on 
direct experimental data from Refs. [SI E] for the dominating resonant contribution, and 
it seems prudent to call for a new R-matrix fit, which is beyond the scope of the present 
work. 

Previous one-zone nucleosynthesis calculations of novae [30] have shown that a 
factor two lower 15 N(p,7) 16 rate results in up to 22% reduction in the final 16 yield, 
depending on the nova temperature. Further implications of the changed 15 N(p,7) 16 
rate are yet to be studied. 

4. Summary 

The 15 N(p,7) 16 cross section has been measured at energies corresponding to hydrogen 
burning in novae. The present data are more precise than previous direct experiments 
[SI E]. They are about a factor two lower than the values adopted in reaction rate 
compilations [7J [8] . 
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Table 1. Spectrum integration and background subtraction. The raw counts in the ROI (9.7-13.5 MeV) and in the background monitoring 
region (14.4-18.0 MeV) are given. For the ratio R^£f°Q , the experimental data are from runs with helium gas in the target. The simulations 
A and B and the adopted uncertainty are explained in the text. The boron background in the ROI (column 8) is obtained by multiplying 
columns 3 and 7. The net counts in the peak (column 9) are obtained by subtracting column 8 from column 2. 



Raw counts Rf£j™ B oron background Net counts 



[keV] 9.7-13.5 14.4-18.0 Experiment Sim. A Sim. B adopted 9.7-13.5 9.7-13.5 A stat A Bo ron 



101 


293 


69 


1.9±0.6 


0. 


70 


1.13 


1.9±1.0 


130±70 


164 


10% 


43% 


122 


355 


31 


3.3±1.0 


0. 


75 


1.37 


3.3±1.0 


100±40 


252 


7% 


14% 


131 


662 


79 


1.8±0.5 


0. 


77 


1.33 


1.8±1.0 


140±80 


522 


5% 


15% 


141 


1703 


172 


1.1±0.2 


0. 


80 


1.44 


l.lil.O 


190±170 


1510 


3% 


11% 


151 


2739 


433 


1.6±0.2 


0. 


80 


1.46 


1.6±1.0 


700±400 


2047 


3% 


21% 


188 


12126 


2895 


1.4±0.3 


0. 


81 


1.46 


1.4±1.0 


4100±2900 


8064 


1% 


36% 


201 


1300 


222 


1.40±0.05 


0. 


84 


1.52 


1.4±1.0 


310±220 


990 


4% 


22% 


210 


32569 


3836 


1.42±0.03 


0. 


80 


1.49 


1.4±1.0 


5400±3800 


27120 


1% 


14% 


216 






1.45±0.09. 


0. 


80 


1.45 












221 


6360 


902 


1.72±0.07 


0. 


78 


1.46 


1.7±1.0 


1600±900 


4805 


2% 


19% 


229 


1930 


98 


2.9±0.2 


0. 


80 


1.43 


2.9±1.0 


280±100 


1649 


3% 


6% 


238 


1517 


33 


2.2±0.8 


0. 


77 


1.51 


2.2±1.0 


70±40 


1443 


3% 


2% 


250 


958 


18 


6±4 


0. 


81 


1.38 


6±4 


110±70 


847 


4% 


8% 



At E p = 216 keV, argon gas has been used instead of helium. 
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Table 2. Systematic uncertainties and their effect on the S-factor data. 



Source of the uncertainty 


Details found in 


Effect on S-factor 


Target density 

O ------- -~- j 


Ref. 125] 

I i 


3.2% 


15 N isotopic ratio 


Refs. [2S1 122] 

E ' 1 i 


3.0% 


Beam intensity 


Refs. [23 [23 

I V I i 


1.0% 


Effective energy 


Ref. [23] 


0.7% - 3.2% 


7-ray detection efficiency 


Sec. |2.2| 


3.0% 


'Y-r&v canture to excited states 


Sec 12 21 


1.9% 


7-ray angular distribution 


Sec. |22| 


4.0% 


11 B(p,7) 12 C background 


Sec. El 


1.8% -43% 


Total systematic uncertainty: 




8% - 44% 



Table 3. Effective center-of-mass interaction energy E , S-factor data, and relative 
uncertainties. The systematic uncertainty due to the boron background subtraction 
has been derived in tableland is repeated here (column 5). The boron uncertainty is 
already included in the total systematic uncertainty given below (column 4). 

[keV] [keV barn] statistical total systematic systematic (boron) 



90.0 


38.4 


14% 


44% 


43% 


109.3 


44.4 


11% 


16% 


14% 


118.5 


47.0 


6% 


17% 


15% 


127.9 


55.4 


3% 


13% 


11% 


136.6 


57.6 


4% 


22% 


21% 


173.0 


72.2 


2% 


37% 


36% 


183.2 


86.1 


4% 


24% 


22% 


192.3 


83.8 


1% 


16% 


14% 


202.8 


85.9 


2% 


20% 


19% 


210.3 


99.9 


3% 


9% 


6% 


219.4 


110.4 


3% 


7% 


2% 


230.0 


120.9 


5% 


11% 


8% 



